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Aims. Recent observations of the line of sight velocity profile of elliptical galaxies have furnished controversial results with some 
works favouring the presence of a large amount of dark matter in the outer regions and others arguing in favour of no dark 
I ■ matter at all. In order to shed new light on this controversy, we propose here a new phenomenological description of the total 
mass profile of galaxies. 

Methods. Under the hypothesis of spherical symmetry, we assume a double power - law expression for the global mass - to light 
(hereafter M/L) ratio T(r) = M(r)/L{r) = To(r/ro)"(l + r/ro)'' with (To,ro) scaling parameters and (a,/3) determining the 
inner and outer slope of the M/L ratio. In particular, T oc for r/ro << 1 so that a = Q mimics a constant M/L ratio in 
the inner regions, while, for (r/ro >> 1), T cx r"^'' thus showing that models with a + /3 = have an asymptotically constant 
Q M/L ratio. A wide range of possibilities is obtained by varying the slope parameters (a, (5) in the range we determine on the 
^ \ basis of physical considerations. Choosing a general expression for the luminosity density profile j(r), we work out an effective 
CZ3 . galaxy model that accounts for all the phenomenology observed in real elliptical galaxies. We derive the main dynamics and 
gravitational lensing properties of such an effective model. 

Results. We analyze a general class of models, able to take into account different dynamical trends. We are able to obtain 
, . analytical expressions for the main dynamical and lensing quantities. We show that constraining the values of a + (3 makes it 
I^nJ ' possible to analyze the problem of the dark matter in elliptical galaxies. Indeed, positive values of a + /3 would be a strong 
evidence for dark matter. 

Conclusions. Finally we indicate possible future approaches in order to face the observational data, in particular using velocity 
dispersion profiles and lensed quasar events. 
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1. Introduction 

Applying the virial theorem to the Coma cluster to esti- 
mate its mass, Zwicky (1937) provided the first evidence 
of the presence of significant amounts of dark matter. 
Notwithstanding the great theoretical, computational and 
observational efforts spent by both astronomers and par- 
ticle physicists, the nature, the properties and the dis- 
tribution of the elusive component remain still unsolved 
problems more than 70 years after the pioneering work of 
Zwicky. Solving these problems stands as one of the most 
challenging yet fascinating issues of modern astrophysics. 

Although detectable only through its gravitational ef- 
fect, dark matter leaves its imprint on a wide variety of 
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observational probes. Rotation curves of spiral galaxies 
(see, e.g., Sofue & Rubin 2001 for a comprehensive re- 
view) are a textbook example. The presence of HI and Ha 
gas makes it possible to reliably determine the kinemat- 
ics and dynamics of these systems. Modelling the galaxy 
with the luminous and the gas components only is unable 
to fit the observed rotation curve unless the mass - to - 
light (hereafter M/L) ratio is unacceptably higher than 
what is predicted on the basis of stellar population syn- 
thesis models. Going up to the cosmological scales, the 
remarkable success of the concordance ACDM model (ac- 
cording to which the cosmological constant A and the cold 
dark matter are the main contributors to the total energy 
density budget) in fitting almost the full set of available 
astrophysical data ( Spergel et al. 2003 IRiess et al. 2004 |) 
provide a further strong evidence in favor of dark mat- 
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ter. Moreover, in this same framework, the presence of 
CDM is mandatory in order to reproduce the observed 
matter power spectrum inferred from large galaxy surveys 
Pawkins et al. 20031 [Pope et al. 2004| . 

However, despite the great success on cosmological 
scales, the concordance model is affected by serious 
shortcomings on the galaxy scales. Indeed, numerical 
simulations of structure formation in the ACDM scenario 
furnish some hints on the mass profiles of dark matter 
haloes. While there is a general consensus that relaxed 
haloes exhibit a density profile that is well described by a 
double power law with outer asymptotic slope —3, there 
is still controversy about the value of the inner asymp- 
totic slope with proposed values mainly in the range 
~ 1.0 — 1.5, predicting the so called cuspy mass density 
models (see INavarro et al. 19961 INavarro et al. 19971 
IMoore et al. 1998|. [G hT gna et al. 2000} IPower et al. 2003| 
INavarro et al. 2004|) . The first results indicated 
an inner slope equal to 1 (|Navarro et al. 1996} 
INavarro et al. 1997| . on the contrary other studies gener- 
ated a more steeper slope equal to 1.5 ([Moore et al. 19981 



Ghigna et al. 2000 



Fukushige fc Making 2001 ) 



Whichever is the actual value of the inner slope, 
observations of rotation curves strongly disfavor the 
presence of any cusp giving, on the contrary, strong 
support to shallower density profiles (the so called cored 
models with an inner slope < 1) in fierce disagreement 
with results of simulations ( McGaugh &: de Blok 19981 
[Binney fc Evans 2001[ Ide Blol» 



c et al. 2001]) . It is worth 



stressing that the usual modeling of spiral galaxies does 
not take into account the presence of spiral arms and 
the possible clumpy distribution of the gas component. 
Both these effects may introduce non circular motions 
causing systematical errors in the determination of 
the rotation curve or in inferring constraints on the 
mass model from the fit to the data. Although recent 
higher resolution simulations seem to partially alle- 
viate both the cusp and the substructure problems 
([Maccio et al. 2006al IMaccio et al. 2006b[) . several alter- 
native scenarios have been proposed in order to solve 
these issues ranging from modifications of the fundamen- 
tal properties of dark matter particles to changes in the 
laws of gravity ([Carroll et al. 2005[) . A different kind of 
approach is followed by Miicket & Hoeft (2003), that 
solving the Jean's equation with suitable assumptions on 
the gravitational potential found an inner slope < 0.5. If 
confirmed, these contrasting results may offer the chance 
to falsify the CDM paradigm on galaxy scales or, at 
least, can furnish a fundamental support to drive N-body 
simulation analysis. 

It is therefore interesting to look at elliptical galaxies 
in order to avoid some of the observational problems de- 
scribed above and to investigate whether dark matter is 
indeed ubiquitous as expected. Unfortunately, here the sit- 
uation is complicated by both observational and theoreti- 
cal difficulties. On one hand, the lack of a well understood 
mass tracer makes it difficult to measure kinematics out to 
the largest radii where dark matter is supposed to domi- 



nate. On the other hand, interpreting the results is compli- 
cated by the strong degeneracy between mass model and 
anisotropy in the velocity space. As a consequence, the 
presence of dark matter in elliptical galaxies is still a mat- 
ter of debate. Using different tracers and modeling tech- 
niques. Mould et al. (1990) and Franx et al. (1994) found 
positive evidences of massive dark matter halos around 
early type galaxies, while no dark matter is needed to fit 
the measured velocity dispersion in some of the galax- 
ies considered by Bertin et al. (1994) and Gerhard et al. 
(2001). Satellite dynamics ( [Romanowsky et al. 2001[ ) and 



X - ray mass estimates ([Loewenstein fc White 19991) adds 
further support in favour of dark matter, while opposite 
results have been recently obtained using planetary nebu- 
lae to probe the outer regions ( [Romanowsky et al. 2003| . 
Furthermore, Mamon & Lokas (2005a, b) have recently 
shown that NFW - like density profiles are not able to re- 
produce the global mass profile of elliptical galaxies whose 
internal kinematics turn out to be well fitted by constant 
M/L ratio models thus arguing against the presence of 
dark matter haloes. 

In order to improve our knowledge of dark matter 
haloes, one has to look for alternative mass probes such 
as gravitational lensing. Born as mere scientific curiosity 
predicted by general relativity, gravitational lensing has 
nowadays given rise to a full sector of modern astronomy 
given its versatility that renders it an ideal tool to inves- 
tigate a wide range of astrophysical phenomena (see, e.g., 
[Schneider, Ehlers fc Falco 1992[ and IPetters et al. 20011 
for illuminating textbooks). Since elliptical galaxies rep- 
resent ~ 80% of lens galaxies in multiply imaged quasar 



systems ( [Fukugita fc Turner 1991 ), reconstructing the 
mass model for these systems gives interesting constraints 
on the dark matter haloes surrounding the lens (see, 
e.g., IKochan ek 1995, Romanows ky fc Kockanek 1997[ 



IKeeton et al. 1998. jTreu fc Koopmans 2004 for a far to 
be exhaustive list of references). Galaxy - galaxy lensing 
([Fisher et al. 2000[ [Guzik fc Seljak 2002[ ) is another 



opportunity to investigate dark matter resorting to the 
weak lensing regime. 

While the use of lensing data alleviates the problem 
with observational issue, there is still a theoretical prej- 
udice concerning the choice of dark matter profile, which 
can crucially infiuence the results. Even if there is a long 
and illustrious tradition of mass models for early type 
galaxies, ranging from distribution function based ones, 
to mass density profile or gravitational potential poten- 
tial based ones, we present a new parametrization of the 
mass profile, describing the M/L ratios as a double power 
law. Thus, we propose here to replace the actual elliptical 
galaxy with an effective one with a mass profile M{r) given 
as T(r)xL(r), with L{r) the deprojected luminosity distri- 
bution. A key role is played by the global M/L ratio T(r) 
which we parameterize as T(r) = To(r/ro)"(l -I- r/ro)^ 
with To a strength parameter, tq a length scale, while a 
and (3 determines the asymptotic behaviors of the M/L 
ratio for r/r^ << 1 and r/rp >> 1, respectively. The 
ansatz adopted for T(r) has the virtue of being quite sim- 
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pie so that many dynamical and lensing quantities may 
be expressed analytically or in terms of special functions. 
Nevertheless, it is also quite versatile giving the possibil- 
ity to mimic the main features of a wide range of modelfl 
For instance, constant M/L models are obtained by set- 
ting (q!,/3) — (0,0), while models with a + fi > Q mimic 
the effect of dark matter haloes on the global M/L ratio. 
Moreover, in the inner regions the main kinematic and 
dynamical quantities does not depend on /3, and if we set 
a = will only depend on the M/L strength Tq and 
the parameters which characterize the luminosity profile 
in agreement with the tentative conclusion of Mamon & 
Lokas (jMamon fc Lokas 2005aP that mass follows light in 
the interiors of elliptical galaxies. 

The paper is organized as follows. In Sect. 2, we in- 
troduce our parametrization of the global M/L ratio dis- 
cussing the role played by the slope parameters (a, (3) and 
the asymptotic behaviours. The basic properties of the 
resulting effective galaxy model are presented in Sect. 3 
where we also compare our model with some widely used 
alternative descriptions of elliptical galaxies. Dynamical 
and lensing properties of the model are presented with 
great detail in Sects. 4 and 5 respectively. A summary of 
the main results and some perspectives on constraining 
the model parameters through observations are presented 
in the concluding Sect. 6. 

2. A variable M/L ratio 

The approach to the problem of dark matter in ellipti- 
cal galaxies is a two step process, where, first, the pro- 
jected light distribution is deconvolved to a three dimen- 
sional mass density by assuming a constant M/L ratio 
and second, a dark halo mass profile is added to the lu- 
minous component. The theoretically predicted luminosity 
weighted velocity dispersion (projected along the line of 
sight) is then compared to the observational data, in order 
to determine the stellar M/L ratio and the parameters of 
the halo model. A successful result of such a fitting proce- 
dure is considered as an evidence of the presence of dark 
matter, while the possibility of fitting the same data with a 
constant M/L model and no dark halo argues against the 
presence of this unseen component. Needless to say, the 
procedure crucially depends on the choice of the model 
for the dark matter distribution. 

In order to alleviate such a problem we work out a fam- 
ily of galaxy models which is able to reproduce, through an 
appropriate parametrization, a wide variety of behaviours 
in the galaxy mass density distribution, including both 
the situation where the dark matter is dominating as well 
the situation where the stellar component is predominant 
As a consequence such a kind of models can be used to 
analyze and shape a whole sample of galaxies in a ho- 
mogeneous way. To this end, it is worth noting that in 

^ Another possible phenomenological approach has been re- 
cently investigated by ( |Napolitano et al. 2004[ ), but they only 
considered an empirical gradient of the M/L ratio assuming 
an a priori model for the dark matter halo density profile. 



constant M/L models, the mass follows the light, while 
this does not hold anymore in presence of a dark halo. 
Actually, defining T(r) = M{r)/L{r) with M{r) and Hr) 
the total mass and luminosity respectively within a dis- 
tance r from the galaxy center, it is straightforward to 
understand that, in presence of a dark matter halo, M{r) 
still increases for values of r beyond the visible edges of 
the galaxy so that T(r) turns out to be monotonically in- 
creasing function in the outer regions. On the opposite, 
in absence of a dark matter halo, M{r) follows L{r) thus 
leading to T(r) ~ constant. Note, however, that, strictly 
speaking, having a constant M/L is not a sufficient probe 
against the presence of dark matter. Actually, it is possi- 
ble that T(r) is a weak function of r over the range probed 
by the data so that detecting its variation is quite hard 
if the data are not of very good quality. Should, however, 
T(r) ~ To ~ with T^, the stellar M/L ratio predicted 
by population synthesis models, then one could indeed ar- 
gue that there is no need of dark matter. As a consequence, 
we may reconsider the problem of dark matter in elliptical 
galaxies as that of directly understanding whether T(r) is 
constant or not and, if yes, what is its trend as a function 
of the radius. Thus, we assume a phenomenological ansatz 
for the M/L ratio that is able to smoothly interpolate be- 
tween the two opposite cases. The only other ingredient 
needed to model the galaxy is then the luminosity density. 
Although it is in principle possible that the stellar M/L 
ratio changes with r because of a distance dependent dis- 
tribution of massive stars, it is worth stressing that color 
gradients are often too low to motivate a strong varia- 
tion of T^. As a consequence, we may safely consider an 
increasing M/L as a clear signal of a dark matter halo. 

Assuming spherical symmetry, we propose the follow- 
ing expression for the M/L ratio : 

with To a scaling M/L ratio, ro a reference radius and 
(a, 13) two slope parameters that we will constrain later 
on the basis of physical considerations only. Note that for 
(a,/3) — (0,0) the model reduces to the constant M/L 
case with T(r) = To. A no dark matter model is then 
obtained setting (a,/3) = (0,0) and To — T^,. 

The slope parameters (a, (3) drive the asymptotic be- 
haviors of T(r) : 

{j.a r/ro << 1 

(2) 

so that it is clear that a determines the slope of the M/L 
ratio in the inner regions, while (3 enters the determina- 
tion of the outer slope. As shown in Fig. 1, the reference 
radius ro marks the transition from the inner to the outer 
asymptotic slopes. 

Eq. (12) helps drawing some useful qualitative consid- 
erations to constrain the parameters {a, (3). First, let us 
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Fig. 1. Global M/L ratio T as a function oi y = r/r^ 
setting To = 4, /3 = 0.3. The inner slope a takes values 
from —0.3 to 0.3 (in steps of 0.1) from the lightest to the 
boldest curve. 



remind that, since both M{r) and L{r) vanish at r — 0, 
we expect that T(r) could diverge, take a constant value 
or converge to 0. In Fig. [T] we plot T as a function of the 
radial coordinate r, fixing (3 and allowing a to vary and 
assume both negative and positive values. An inspection 
of the figure shows us that negative values of a give diverg- 
ing values of T for r — s- 0, on the contrary null or positive 
a give a T that becomes finite or null. Having T(r) — s- cx) 
as r ^ could seem unrealistic since one could argue for a 
huge content of dark matter in the inner regions, but this 
is not the case, since it is possible that L{r) more 
rapidly than M(r). 

On the other hand, it is easy to see that models with 
a + (3 — describes galaxies with finite total mass since, in 
this case, both AI{r) and L(r) for r >> tq are constants. 
Moreover, since M(r) is asymptotically increasing or fiat, 
while L(r) is constant, T(r) should be an increasing or flat 
function of r for r » tq thus leading to the constraint 
a + P > 0. As we will see later, we will recover these same 
constraints from a different and more physically motivated 
perspective. 

As a final remark, let us remember that we have 
assumed spherical symmetry. Actually, real galaxies are 
hardly spherically symmetric, but are intrinsically ellip- 
soidal or triaxial systems. However, an ellipsoidal or tri- 
axial version of the model we will investigate may be 
obtained by replacing in Eq. ^ the spherical radius 



3. Basic properties 

The ansatz ([T]) is the main ingredient of our approach 
to the dark matter problem in elliptical galaxies. In or- 
der to complete the description of a galaxy, we need 
to choose a luminosity density profile that, when pro- 
jected on the plane of the sky, is able to fit the ob- 
served surface brightness. The ideal route to follow should 
be to deproject the Sersic profile (jSersic 1968.) since this 
is known to well fit the surface brightness of ellipti- 
cal galaxies (jCaon et al. 1993|. IGraham fc Colless 1997|. 
[Prugniel Sz Simien 1997 ). Although this is mathemati- 
cally possible ( Mazure fc Capelato 2002[ ), the result is ex- 
pressed in terms of unusual special functions and is there- 
fore of little practical utility. A different way is to choose a 
luminosity density profile whose two dimensional projec- 
tion mimics the de Vaucouleurs law llde Vaucouleurs 1948[) 
since this profile also gives a good fit to elliptical galax- 
ies. A versatile family of models suitable for our study is 
provided by the Dehnen (1993) models whose luminosity 
density is : 



(3) 



with Jo and a scaling quantities and 7 a positive slope pa- 
rameter. Note that the model reduces to the Jaffe (1983) 
and Hernquist (1990) ones for 7 = 2 and 7 = 1 respec- 
tively. The luminosity density scales as r^"' at the cen- 
tre, while asymptotically decreases as r~^. The constraint 
7 < 3 ensures that the model is physically well behaved, 
but, actually, only models with 7 < 2 are interesting 
since their surface density fits well the de Vaucouleurs 
law provided the scale radius a is expressed in terms of 
the effective radius Re. In particular, for 7 = 1, it is 
a = 1.33(1 -I- V2)~^Re, while similar relations hold for 
other values of 7 (jPehnen 1993) . Models with 2 < 7 < 3 
are still meaningful, but have some disturbing features 
such as an infinite central potential and a diverging cen- 
tral velocity dispersion. In the following, we will set 7 as 
a free parameter, but in the figures we will set 7 = 1 in 
order to recover the widely used Hernquist model. 

Being the model spherically symmetric, it is immediate 
to evaluate the luminosity density. Inserting Eq. ^ into 
the standard relation 



(x2 + J/2 + ^2/^2)1/2 ^i^g L{r) = / 4TTr'^j{r')dr' 



r with the elliptical m = {x'^ + y'^ + 
triaxial one m — {x^ + y^/p^ + jq^Y^"^ with q and 
•p the major- to- minor and intermediate - to - minor ax- 
ial ratios ( [Binney fc Tremaine 1987| . Models which result 
from equation ([1]), after such replacement, are similar to 
the spherical ones in some important properties such as 
mass profiles. Moreover, it is worth remembering the exis- 
tence of a degeneracy among mass, anisotropy and fiatten- 
ing that makes it possible, at least in principle, to mimic 
a flattening of the system with an anisotropic velocity dis- 
tribution ( [Magorrian fc Ballantyne 200 Ij ). Finally, using 
spherical symmetry allows us to express the main quanti- 
ties of interest in terms of analytical or special functions. 



we easily get : 

i(r) =it„ta;3-T(l+x)-^-3 

with X = r I a and Ltot the total luminosity given by : 
47rjoa^ 



L 



■ 7 



(4) 



(5) 



Note that the constraint 7 < 3 ensures that the total 
luminosity is flnite and that L[r) does not diverge at the 
centre as it is physically reasonable. 
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3.1. Mass profile 

Having chosen a model for the luminosity density, the 
mass profile of the galaxy is immediately obtained from 
the definition of M/L ratio and the ansatz ^ : 



M{x) 



To Ltot (1 + f 1 + — 



(6) 



with xq = ro/a. In particular, we note that, the larger 
are /3 or a, the greater is the value of the total mass. For 
a + (3 = we have a finite mass. 

Eq. ([6]) makes it possible to put physically motivated 
constraints on the slope parameters {a, P). First, the mass 
profile must vanish at the centre so that we get : 

M(r = 0) = 0^3 + a- 7>0 4=> a > 7 - 3 . (7) 

Let us now consider the asymptotic behaviour at infinitjQ : 



00 for a + /3 > 

lim M (r) <( Tq Ltot for a + /3 = 
for a + /3 < 



(8) 



Needless to say, a null total mass is unphysical so that 
Eq. ^ leads to the constraint : 



a + (3>0 



a>-(3. 



(9) 



Note that the case a + (3 — Q individuates models with 
finite total mass as we have argued above directly from 
Eq. ©. 

3.2. Mass density profile 

For spherically symmetric models, the mass density may 
be evaluated differentiating the mass profile. Inserting 
Eq.® into the usual relation 



p(r) 



1 dM 1 dM 



Anr^ dr Ana^x'^ dx 
we easily get : 



p(r) 



_ Tq Ltot a-7 



/3-1 



First, we note that models with a cusp are obtained for 
a — 7 < 0, while a flat inner core is present for models 
with a — 7 = 0. Moreover, since the mass density cannot 
increase with r, we get the constraint : 



a — 7<0— !-a<7 



(13) 



Let us study now the behaviour in the opposite limit, 
i.e. for a: >> 1. In this case, we have to consider separately 
models with infinite or finite total mass. In the former 
case, a + f3 ^ and we get : 



p{x) cx x"+'^"3 for X » 1 , 
while for models with a + /3 = 0, it is : 
c^'' for 3 — 77^ axo 



p{x) cx 



n-5 



for 3 — 7 = axQ 



(14) 



(15) 



Imposing that /? is a decreasing function at infinity makes 
it possible to further constrain the slope parameters. 
Considering Eqs. (fT4|) . we thus get : 



(16) 



a + /3 - 3 < ^ a < 3 - /3 



3.3. Surface mass density 

It is worth deriving the surface mass density of the model 
since this will enter the derivation of the lensing quantities 
which Sect. 5 is devoted to. To this aim, we start by the 
usual definition : 



S = 



P{ri,r2,r3)dr3, 



(17) 



with (ri,r2) cartesian coordinates in the plane orthog- 
onal to the line of sight and along the line of sight. 
Inserting the mass density (IT5|) in Eq. pT|) and changing 
to cylindrical coordinates, one gets (see Appendix A in 
ISand et al. 20041) : 



m = 



J^tot i / 



2 7rr, 



;^l + Q — 7 



csc(6i) 



2+Q-7 



[l+xy-^ 1 + - ^(^)(10) X (1 + C csc(0))-4+7 (a-Q + ^ csc(0))-i+^ 



with 

V = (a+/3)a;2 + [3 + a{\ + xq) + /? - 7] .T+(3+a-7)a:o -(U) 

The mass density in Eq. PH)) must be physically reasonable 
so that we can further constrain (a, [3). To this aim, let us 
consider its asymptotic behaviours. In the inner regions, 
it is : 



X [xq (3 + a — 7) + (3 + a + a + /3 — 7) ^ csc(^?) 
+ {a + P) csc{ef]de. (18) 



p{x) (X x" for X «1 . (12) 

^ Although a finite total mass is a desirable feature, there 
are plenty of interesting models with a formally infinite total 
mass that are typically truncated by hand at the virial radius 
thus obtaining a finite total mass. 



where ^ = R/a, with R — \Jr\+r\. Eq. lfTS]) may be 
analytically evaluated only in the case xq — \ otherwise 
numerical integration is needed. However, also for = 1, 
the result turns out to be a rather messy combination 
of special functions so that it is preferable to resort to 
numerical integration also in this case. 

The asymptotic trends of the surface mass density de- 
pend on which is the range a lies in. If 7 — 2 < a < 7 — 



See Eq. (|32|l to motivate the lower bound, more stringent 
than Eq. ((Tjl 
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— > ^" for ^ — > 0, while, S asymptotes to a con- 
stant value in the same limit for (7 — 1) < a < 7. These 
regimes were qualitatively analyzed in the previous sec- 
tion and can be inspected in Fig. [U where it is clear how 
the inner trend of the M/L ratio is crucially dependent 
on the choice of a. On the other hand, S oc ^"+/3-2 

^ 00, which is expected given the outer slope of the 
three dimensional mass density. 

Further constraints on the slope parameters (a, (3) may 
be obtained considering the above asymptotic scaling of 
E(^). Indeed, imposing that E(^) is always decreasing, 
(a, (3) must satisfy the constraints a — 7 -|- 1 < and 
a + (3 — 2 < 0. In particular, for 7=1, the first constraint 
translates to a < 0, and an asymptotic flat S(^) for ^ ^ 
is obtained when is a > 0. Thus, a crucially influences 
the inner trend of E similarly to what done for T. The 
other constraint further reduces the range of variability of 
a and /3 in Eq. (|16p , giving the condition a < 2 — (3. 



3.4. A general comment on our effective model 

A general comment is in order here. The philosophy moti- 
vating our approach to the dark matter in elliptical galax- 
ies is favored by the phenomenology of the problem. We 
have adopted an ansatz for the M/L ratio and a model 
for the luminosity distribution and are now working out 
the basic properties of a system that may be described by 
these two assumptions. Strictly speaking, this system does 
not represent the actual galaxy, but it is an effective model 
of the galaxy. Should dark matter be indeed present, the 
galaxy must be correctly modelled as a two components 
system with a dark halo embedding the luminous compo- 
nent. However, from the point of view of an observer, the 
dynamical quantities that can be measured and compared 
with the theoretical predictions (such as the rotation curve 
and the velocity dispersion) are determined by the total 
mass and density profile only. As a consequence, describ- 
ing the galaxy with an effective model is phenomenologi- 
cally equivalent to modeling the system with two different 
components. Actually, one could fit an effective galaxy 
model and then decompose it in a luminous component 
plus a dark matter halo determining the density profile of 
the halo as ph — p — Ti,j{r) with p given by Eq. (fTO|) . In 
the phenomenological approach we are pursuing here, we 
only need to study the dynamical and lensing properties 
of the effective model defined by Eqs.(IT]) and ([3]), while 
disentangling the two components is not necessary. 

It is, nevertheless, interesting to compare the relative 
contributions of ph and = T^,j(r) to the total mass 
density p. Depending on the values of the parameters 
{a,f3,Xo), the two components may furnish comparable 
contributes to p for medium radii (x around 1 or few 
units), while, as expected, it is the dark matter halo to 
dominate the outer mass budget. Nevertheless, given the 
way we have assigned the model, the dark matter density 
has the same asymptotic behaviours of the total one, i.e. 
Ph oc x"^'' for X ^ and ph oc for x 00, thus 



the inner slope of the dark matter density is determined 
by the luminous component through the value of 7. As a 
consequence, the constraints on the model parameters a 
and 7 provide information about the competition between 
the luminous and dark matter of galaxies, especially in the 
inner regions where the baryons are supposed to be dom- 
inating (jMamon fc Lokas 2005"a|) . 

3.5. Local M/L ratio 

A general comment is in order here. Eq. ([T]) defines the 
M/L ratio as M{r)/L(r). This is also referred to in liter- 
ature as global M/L ratio, while the quantity Tioc{r) = 
p{r)/j{r), with p{r) and j(r) the mass and luminosity den- 
sity, is dubbed local M/L ratio. 

Now it is interesting to compute the local M/L ratio 

as : 



Tioc(a;) 



~ 3{x) (3 - 7)xo(l + x/a;o) 



T{x)V{x) 



(19) 



Apart the constant numerical factor (that is positive defi- 
nite for 7 < 3), the local M/L ratio is similar to the global 
M/L ratio, but is no more a double power law because 
of the multiplicative term V{x)/{1 + x/xq). Note that, 
in the inner regions (x << 1), both T(a;) and Tioc{x) 
scale as x". On the other hand, the behaviour at infin- 
ity is radically different. To this end, let us note that, 
for a; >> 1, T^oc scales as 0;"+''+^, whereas Eq.Q shows 
that T oc x"~^^ . As a consequence, models with finite to- 
tal mass (i.e. a + P ~ 0) have an asymptotically constant 
global M/L ratio, while a divergent local M/L ratio. We 
thus argue in favour of using T rather Tioc in order to 
avoid any disturbing divergence. 

3.6. Comparison with a generalized NFW profile 

It is interesting to compare our phenomenological ansatz 
for the M/L ratio with that coming out from an assumed 
dark halo model. To this aim, describing in both cases 
the luminosity component with a Dehnen model, we have 
to choice a dark matter density profile among the various 
proposals available in literature. Motivated by the results 
of a wide set of numerical N-body simulations, it is cus- 
tomary to model dark halos through the so called NFW 
model whose mass density profile reads (Navarro et al. 
1996,1997): 



PNFwir) ^ Ps[ — 
\rh 



-1 / \ -2 

r 

1 + — 

rh. 



(20) 



where and ps scaling values for the radius and the 
density. Following the common practice, it is more con- 
venient to parametrize the model by the virial mass M^ 
and the concentration c = r^/r^ with r„ the virial ra- 
dius. As it is easily seen, pnfw asymptotically scales 
as and in the inner and outer regions respec- 
tively. While there is a wide consensus on the value of 
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the outer slope, a fierce debate is still open on the in- 
ner scaling with values as high as 1.5 (jMoore et al. 19981 
[Ghigna et al. 2000HFukushige fc Makino 200ip . As a con- 
sequence, it is more instructive to consider the generalized 
NFW (gNFW) model defined as : 



PgNFwir) ^ Ps[ — 



-5 / ^ X 5-3 



(21) 



which reduces to the NFW model for 6=1, while 6 = 1.5 
reproduces the results obtained by Moore et al. (1998). 
Note that the total mass is formally infinite, but, as usual, 
we truncate the model at the virial radius and take My as 
finite total mass. 

In the rest of this section, we will compare the Ad/L 
ratio and the mass density profile of our effective galaxy 
with the same properties of a real galaxy modelled with a 
Dehnen-like luminous component and a gNFW halo. For 
illustrative purposes, we also consider the case of a sin- 
gle gNFW component although this leads to less reliable 
results in the inner baryon dominated regions. 

3.6.1. M/L ratio comparison 

The global M/L ratio for the gNFW + Dehnen (here- 
after gNFW-|-D) models may be easily computed as 
'^gNFW+D = [MgNFW + ^ *Lir)]/ L{r) , while for the 
gNFW only model it is simply Tgjypw — MgNFW /L{r). 

Depending on the values of (7,(5), we can have differ- 
ent variation of T g^pw+D with r. For instance, if 5 > 7, 
^ gNFW+D presents a mimimum and increases outwards, 
while in the inner region may eventually diverge. Similar 
trends may be observed also for Tg^pw which is not sur- 
prising given that TgNFw+D{r) -I- Tg^pwir). Our 
phenomenological ansatz for T(r) in Eq.([T]) presents sim- 
ilar characteristics with a playing the role of 5. In partic- 
ular, models with a < presents a minimum and may be 
divergent at the centre. 

The similar role played by a and d may be well under- 
stood considering the asymptotic trends of Tgjv fw+ d and 
^gNFW- First, when 7 < (5, for r ^ both T gpfpw+D and 
'^gNFW scales as r''~^ so that both quantities diverge. On 
the other hand, if 7 > 5, we get T g^pw+Di^) ~ while 
^ gNFw{0) = 0. In the opposite limit, for r ^ 00, we get 
^gNFW+D = ^gNFW = Mfot/Ltot, i-B. an asymptotically 
constant value. 

Our ansatz for T scales as r" in the inner regions and 
as r"~^^ in the outer ones, being therefore able to recover 
the trends of gNFW-|-D and gNFW models. For instance, 
in the case of gNFW-f D model, the inner slope is recov- 
ered for a = 7 — 5 or a = depending on the 7 being 
smaller or larger than 6, while an asymptotically constant 
value is obtained for a -|- /3 = 0. Our model can repro- 
duce asymptotically both flat M/L (for a + (3 = 0) and 



* We note that the gNFW model, also being able to repro- 
duce observed flat rotation curve, has a total finite mass and 
therefore an asymptotically flat M/L ratio. However, this is 
not a serious shortcoming such a profile, since only the be- 



increasing M/L (for a + f3 
respect to the gNFW. 



1) giving more possibilities 



3.6.2. Mass density comparison 

For what concerns the mass density, the NFW+Dehnen 
model presents the following trends 

r^^ for r ^ and 5 > 7 

T for T — ^ and S ^ ^ 
PgNFW+D{r) oc { ~ (22) 



for r 



00 



The gNFW model has the same value for the outer slope, 
while PgNFW oc for r — > 0. 

Since the mass density of our effective model scales as 
j-a-'r for r ^ 0, we are able to recover the same inner 
asymptotic trends of the gNFW-|-D and gNFW model by 
suitably choosing a. Considering the case 7 = 1, the choice 
a = gives p cx r^^ as for the NFW model {5 = 1), 
while for a = —0.5 we get the Moore model {S = 1.5). 
Positive values of a do not recover the trend of gNFW-|-D 
model, since in this case p cx r~'^ with e G [0, 1), and 
PgNFW+D OC r"^. In the case of gNFW model, assuming 
7 = 1, the models with d > 1 {< 1) are recovered when a is 
negative (positive). Finally, we note that, once a has been 
set in order to recover the same inner asymptotic trend 
of a given gNFW-|-D model, the outer scaling may 
be recovered by setting /S = —a since our mass density 



scales as 



„a+/3-3 ; 



in these regions (as stated above when 



we analyzed the M/L ratio). 



4. Dynamical quantities 

After having evaluated the main properties of the mass 
model, we determine here some useful dynamical quan- 
tities. In a first step, we derive the gravitational poten- 
tial discussing its main properties. From an observational 
point of view, however, the rotation curve and the velocity 
dispersion are more appealing quantities since they can be 
compared to the observed data in order to investigate the 
viability of our parametrization of the M/L ratio. 

4.1. Gravitational potential 

Although not explicitly needed to evaluate the dy- 
namical quantities of interest, it is nonetheless use- 
ful to compute the gravitational potential for the 
model we are considering. For spherical systems, it is 
dBinney &: Tremaine l"987| : 



$(r) = ^ GM{r) _ p{r')r'dr' 



haviour in a limited and finite region is important to describe 
the observations and on these scales the model reproduces a 
nearly flat rotation curve. 
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Fig. 2. Normalized circular velocity Vg = ^^jfy for different combinations of the parameters {a, (3, xq). Starting on the 
top, in the panel on the left we consider model with finite mass (a + f3 — 0) and on the right models with flat curves 
{a + /3 = 1). In the bottom, on the left we fix a and change /3 = 0.5, 1.0, 1.5 and on the right a ~ 0, (3 — 0.5 and we 
change xq = 0.5, 1, 2. 



that, inserting Eqs.® and (fTU)l gives: 



$(a;) 



X 
Xo 



(l + j;)3-^ 

with p the set of parameters (a, /?, 7, Xq) and 



Xq 



(23) 



for this class of models simplifies to —GM{x)/x which is 
the well known Keplerian potential of a spherical system 
with finite total mass. 



Finally, we study the asymptotic trends. For x — > 00, 
(j) oc x~^~^°'~^^ ^ thus in order to have a plausible gravi- 
tational potential, the models with a + /3 > 1 are ruled 
out, so that we conclude that the only physical models are 
those with 



a + /3-l<0^a<l-/3. 



(28) 



Note that this expression only holds for models with a + 

13 <1 and is a complicated combination of hypergeometric ^ 2 Circular velocity 
and Appell functions. As an example, we report the result 
for xq = 1 which is the simplest one : 



T{x) = 



7 - 3 - a {a + (3)x ' 

Q + P — 1 



For a spherically symmetric mass distribution the circular 
velocity is simply given by 



(25) 



GM{r) 



which, in our case, becomes : 



_a + f3-2 
having defined : 

Ti{x) = 2Fi{2^a~ (3A- P-l,i-c,- p,^\/x) , (26) 
T2{x) = 2Fi{l-a- P-l,2-a^ P,-\/x) , (27) 
with 2i^i(«,5,c,y) the hypergeometric functioii K is interesting to look at the asymptotic trends : 

2,2+a-7 fQj. 2; << 1 



(29) 



a;2+a-7 (1 + a;)7"3 (1 + x/xoY^ .(30) 



vl{x) (X 



X»+0-l foj. 2; >> 1 



(31) 



dGradshteyn fc Ryzhik l980l ). 

The case a + P = 1 merits some discussion since the 
function turns out to be constant. Since the gravita- 
tional potential is defined up to an arbitrary additive con- 
stant, we may choose this constant in such a way that $(x) First, we note that Vc must vanish at the centre so that 



We use here the same notation for the hypergeometric func- 
tion as in the Mathematica package. 



we get the condition : 

2-|-a-7>0-^a>7-2 



(32) 



C. Tortora et al.; Dynamics and lensing in a phenomenological model for ellipticals 



9 



that is more stringent than Eq.([7|). We consider as phys- 
ically motivated those models having an asymptotically 
flat or decreasing rotation curve. Thus, using Eq. pTjl . we 
obtain the constraint in Eq. 

In particular, models with a + f3 ~ presents a keple- 
rian fallofi^ of the rotation curve consistent with the finite 
total mass, while Vc is asymptotically flat when a + /3 = 1 
with : 



Vc{co, a + /3 = I) 



axQ 



(33) 



The qualitative dependence of Vc on the parameters 
{a,l3,xo) may be derived by Fig. [2] where we plot the 
scaled circular velocity Vg = Vc{x)/vc{x — 1) for some 
representative choices of the parameters (a, (3, xq) and set- 
ting 7 = 1. In the top left panel, we present models having 
flnite total mass (i.e., with a -I- /? = 0) thus showing a ke- 
plerian falloff of the rotation curve. The top right panel, 
instead, refers to models with asymptotically flat rotation 
curve that is with a + fi — 1 = Q corresponding to a M/L 
ratio that is still increasing beyond the visible edge of the 
galaxy as in presence of a dark matter halo. In both cases, 
for a given a;, Wc is a decreasing or increasing function of (3 
depending on x being larger or lower than 1 respectively. 
The situation is reversed when considering the dependence 
on a with the higher values of a leading to larger (smaller) 
circular velocities in the regime a; > 1 (a; < 1). The role 
of (3 is further investigated in the bottom left panel where 
we hold a fixed and change (3 showing that this slope 
parameter determines the transition from asymptotically 
decreasing to asymptotically flat rotation curves. In par- 
ticular, models with flat rotation curves rotate faster than 
those with a keplerian falloff. Finally, in the bottom right 
panel, we investigate the dependence on the scaling radius 
Xq finding out that the larger is Xq, the lower is Vc(x). It is 
worth noting that the wide variety of cases is an evidence 
of the extreme versatility of the adopted parametrization 
even if this originates degeneracies among the parameters 
(a,/3,a;o) that may not be disentangled with observations 
of the rotation curve only. 

Combining all the constraints obtained insofar on 
model parameters, we end up with the following range 
for a : 



max (7 — 2, — /3) < a < min (7, 1 — /3) 



(34) 



which we have graphically summarized in Fig.[3]for models 
with 7 = 1. 

4.3. Velocity dispersion 

While the circular velocity is related to the ordered mo- 
tions in the galactic plane, the velocity dispersion takes 
into account disordered motions and thus is the main 
quantity to describe dynamics of non rotating systems like 
elliptical galaxies. Assuming an isotropic velocity disper- 



2 
1.5 
1 

P 0.5 


-0.5 

-1 




0.5 



Fig. 3. Constraints on model parameters a and /3 for 
7 = 1. The selected region is allowed by a theoretical anal- 
ysis on the model and its main related quantities (mass, 
mass density, surface mass density, gravitational potential 
and circular velocity). The parameters a and 7 are crucial 
in this analysis, in particular to determine the low radii 
trend of the surface mass density. In the limit r — > 0, the 
darker region corresponds to diverging ('cuspy') S, while 
the lighter one individuates models with flat ('cored') S. 

sion tensor, the solution of the Jeans equation gives the 
following general formula ( [Binney fc Tremaine 1987| : 

1 r ,.GM{r') 
j{r) 



'-dr' 



(35) 



that, applied to our particular case, reduces to : 



a'ix) = 



with : 



a Xn 



c^-r {1 + xy'-*{l + x/xo)'^-'^'P{x) 



(36) 



^l+2(a-7)(l + ^)27-7(i ^ ^/^p)2/3-lp(^)^^ (37) 



This integral is only defined for a + (3 < 2 which is not a 
problem given that a nd (3 ranges in the interval defined 
by Eq.(l34|). The general analytical expression is again a 
complicated combination of hypergeometric and Appell 
functions so that a numerical integration turns out to be 
easier to handle in the applications. As an example, we 
only report the case with xq = 1 that is : 



„2(q+/3)-5 



2{a + f3-2) [2{a + P)~5] 
{ 2(2-a-/3)(3 + a-7) + 
[ a{2a + 4/3 - 5) + /3(2/3 - 5)] x I2 } 
where we have set : 

Ji EE 2^1(5 -2a -2/3, 7 -2/3 -27, 6 -2a -2/3, 
J2 = 2^1(4- 2a- 2/3,7- 2/3- 27,5- 2a- 2/3, 



(38) 

-l/x),(39) 
-1/x) .(40) 



A wide variety of cases is obtained varying the parameters 
(a,/3, xo). Note, however, that what is indeed observed is 
not the velocity dispersion (t{x), but rather its luminosity 
weighted projection along the line of sight that we will 
compute later. 



10 



C. Tortora et al.: Dynamics and lensing in a phenomenological model for ellipticals 




Fig. 4. A(Ts as a function of (3^, for x — 0.1 (thin 
dashed), x = 0.5 (thick dashed), x = 2 (thin sohd) and 
a; = 10 (thick sohd). The model parameters are set as 
(a, /3, 7,2:0) = (-0.5,1.0,1.0,1.0). 

4.3.1. Effects of anisotropy 

Eq. ((35|) only holds under the hypothesis that the velocity 
dispersion is the same along the three axes of the coordi- 
nate system. In the case of a spherically symmetric system, 
the relation ag = Ucf, — at holds, where tre, (70 and at are 
respectively the polar, azimuthal and tangential velocity 
dispersion and we can define the anisotropy parameter: 



/3.(r)^l-^, 



(41) 



Fig. 5. Aas as a function of fioM- The values of x, the 
line coding and the model parameters are the same as in 
Fig.El 



with roM the anisotropy radius. Isotropic models are re- 
covered in the limit roM — > cxd. In such a case we can 
extend the definition of anisotropy to isotropy ratio pa- 
rameter as Aas — 100 x [as{x, xqm) / as{x, 00) — 1], where 
xoM = tomIo" In Fig. [5] we plot Aa^ as a function of 
[iQM = Although the results are more complicated, 

it is nevertheless clear that also in this case the anisotropy 
introduces large changes in the velocity dispersion. Aas 
converges to for high xqm (low (J-om) and assumes sat- 
urated and very high values for low xqai (high fioAi)- 



with tJr the radial velocity dispersion; isotropy in the ve- 
locity space means f3cr — (i.e., when ar — at). In the 
general case, Ea. (|35)l becomes: 



a^(r 



V[r )j{r ) ^ — dr , 



ri[r)j{r 

where rj{r) is obtained by solving : 
dlogri 



rflog? 



2/3. (r). 



(42) 



(43) 



To quantify the effect of the anisotropy, we first con- 
sider the simplest model /?o-(r) = const and introduce 
the anisotropy to isotropy ratio parameter Aas = 100 x 
[as{x, /3(j)/o's(a;, 0) — 1], defined as a function of scaled ve- 
locity dispersions In Fig.|4] we plot Act^ as a function 
of [3 a for different values of the dimensionless radius x. 
It turns out that a constant anisotropy significantly al- 
ters the velocity dispersion profile, with positive values 
increasing the velocity dispersion at low and high radii. 
The amount of the relative variation depends on x and 
the model parameters (a, /3, 7, xq). 

A constant anisotropy, although simple and widely 
used, is not a realistic assumption. More realistic 
results are obtained considering the Osipkov- Merrit 
parametrization dOsipkov 1979, .Merritt 1985P : 



(44) 



' OM 



<Js {x) = 



^YY, is a scaled version of the velocity dispersion. 



4.3.2. Projected velocity dispersion 

The velocity dispersion defined by Eq. ((35| is not directly 
measurable. Indeed, when comparing to observations, we 
have to take into account some effects. First, the line of 
sight does not coincide with the radial direction so that 
we measure only the projection of the velocity dispersion 
along the line of sight. Moreover, all the galaxy elements 
contribute to the measured a, but the more luminous is 
the element, the larger will be its contribution to the ob- 
served quantity. As such, we have first to define the ve- 
locity dispersion projected along the line of sight and lu- 
minosity weighted. Denoting by aios this quantity, it is 
Pinney fc Tremaine l"987| : 



'los 



I{R) 



' , ,GM(r)Vr^ - , 
J {r) 2 ^'^ 



(45) 



with I{R) the luminosity intensity. 

As a second step, we have to take into account the 
seeing that makes the intrinsic intensity differs from the 
observed one. We therefore define the seeing corrected ve- 
locity dispersion as : 



^ seeinq (-^) 



Jd^R'P{R- R')I{R')al{R') 
Jd^R'P{R-R')I{R') ' 



(46) 



where P{R — R') is the point spread function taking into 
accont both the atmospheric and the instrument seeing. 
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Fig. 6. Normalized deflection angle (left) and lensing potential -08 (right) as function of ^. The model parameters 
(a, a;o) are set as in the inset of each panel, while 7 = 1 is used for the luminosity profile. 



Finally, we have to consider the effect of the spatial bin- 
ning and the finite width of the slit used to make the 
measurements. Therefore, the measured quantity reads : 



bin 



J^dA'Is{R')al,,^g{R') 



(47) 



where A is the area of the slit and Is is the seeing corrected 
intensity. We do not investigate these observational effects 
since they are strictly related to the observational setup. 
We, however, stress the need to carefully take into account 
these corrective terms in order to not introduce systematic 
errors when comparing the model to the data. 

5. Gravitational lensing 

Considered at its early beginning a little more than a 
scientific curiosity, gravitational lensing has now given 
rise to what may be referred to as lensing astronomy. 
Being directly dependent on the mass rather than the 
light distribution, gravitational lensing has turn out to 
be an invaluable tool to investigate the total density 
profile in early type galaxies at intermediate and high 
redshift. The lensing observable quantities (such as the 
position of the images and their flux ratios) in multi- 
ply imaged sources (quasars or galaxies) make it possi- 
ble to constrain the mass profile. Combining lensing re- 
construction with a measurement of the velocity disper- 
sion helps breaking degeneracies inherent in lens modeling 
thus shedding a powerful light on their dark matter con- 
tent dTreu k Koopmans 2002al|Treu fc Koopmans 2004| . 
Understanding the mass distribution of the lens plays a 



key role also in the determination of the Hubble con- 
stant Ho through the time delay method (IRefsdal 1964p . 
In particular, Kochanek (2002) has shown that the pre- 
dicted delay between two images, and hence the derived 
value of Hq , is primarily governed by the average surface 
mass density in the annulus defined by their radial dis- 
tances from the lens center with a small correction taking 
into account the slope of the profile in this annulus. A 
careful determination of mass profile is therefore manda- 
tory to get a reliable estimate of Hq and hence of the 
distance scale. Lens modeling plays a vital role also in 
cosmological applications of gravitational lensing. Since 
early -type galaxies are expected to represent at least 80% 
of lenses (jTurner et al. 19841 [Fukugita fc Turner 199l| ), a 
correct description of their luminous and dark component 
is an essential ingredient in every attempt to constrain cos- 
mological parameters through statistics of lens systems. 

Motivated by these considerations, we complement our 
presentation of our new phenomenological model by a de- 
tailed investigation of its lensing properties. 

5.1. Deflection angle and lensing potential 

For a spherically symmetric model, the deflection angle a 
may be easily computed as : 

2a 

d(0 = — / Kiu)udu = aiXas{C,p) (48) 
s Jo 

where the normalized deflection angle as depends on the 
dimensionless radius ^ — R/a and the model parame- 
ters p = (a, /?, 7, xo), while ai is the value of the de- 
flection angle for ^ = 1. The normalized surface mass 
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density k is defined as the ratio E/Ecrit between the 
surface mass density E and the critical surface density 
T,crit = c^Ds/i^TrGDaDds) with D^, and D^s the 
source, lens and lens -source angular diameter distance^ 

Although quite complicated analytical expression are 
possible for as, we prefer to not report them here and 
discuss the results looking at the left panels in Fig.[6l In 
particular, in the top left panel, we consider models with 
an inner core (a = 0) and find out that the larger is /?, 
the higher (the lower) is the deflection angle for ^ > 1 
(^ < 1). Note that, in the outer regions, as{^) is always 
larger than what is predicted by a constant M/L model. 
In the bottom left panel, we show the effect of changing 
a while holding f3 fixed. It turns out that increasing a 
strongly lowers the deflection angle in the outer regions 
(^ > 1) with only a minor effect for ^ < 1. These trends 
are consistent with the asymtptotic scalings as oc ^°'-f+'^ 
for ^ << 1 and oc for ^ >> 1. Note that, as 

expected, for models with finite total mass (a + /3 = 0), we 
get as (X as for a point mass lens. Finally, we note that 
as is an increasing function of xq which can be explained 
by noting that the larger is xq, the more the model is 
concentrated and hence more mass is within a circle with 
radius ^ and hence the larger is the deflection angle. 

The spherical symmetry of the model makes it pos- 
sible to straightforwardly compute the lensing potential 
integrating the deflection angle thus obtaining : 



(49) 



sin( 







(52) 



V' = a / aiOdC = V'l X (C; p) 



with ips the scaled potential and i/^i — — 1). The 
complete expression of ips is too cumbersome and will be 
not reported here. However, we plot ips as a function of ^ 
for some representative models in the right panels of Fig. [5] 
showing that ips is an increasing function of p). For ^ — > 0, 
Ips oc while Ips oc ^"+'^ as ^ ^ oo. Note that both 

these trends should be derived directly from those of as 
given how as and ips are related. 

The deflection potential is an essential ingredient to 
write down the time delay function and hence estimate 
the time delay between the images in the multiply im- 
aged quasars. Adopting polar coordinates (i?, 9) in the 
lens plane with 6 measured counterclockwise from North, 
it is : 



At 



-ds)+les-H^,o) 



,(50) 



where {^,9) and (Cs,^s) are the images and source po- 
sitions, h the dimcnsionless Hubble constant and rioo = 
{l + Zd)/cxDd.Ds/Dds with Zd the lens redshift. According 
to the Fermat principle, the images form at the extrema 
of the time delay function so that the lens equations are : 



^-6COS(0-0S) = 



aiO 



(51) 



In the rest of the paper we adopt the flat concordance 
ACDM model with {nM,^iA,h) = (0.3,0.7,0.72) consistent 
with the WMAP measurements ( |Spergel et al. 2003[ ). 



Eg. ([5^ reduces to an identity for ^s = 0, i.e. the lens and 
the source are perfectly aligned. In this case, the image 
is the well known Einstein ring with radius obtained by 
solving Eq. ([5T|l for ^s — 0. Alternatively, Eq. lfSS]) is solved 
hy 9 — 9s + mr with n = 0, 1. Inserting these values in 
Ea. (|5ip and solving with respect to ^, we get three images. 
Two of them are disposed symmetrically with respect to 
the lens centre, while the third one lies very close to the 
centre (with ^ ~ 0) and turns out to be unobservable. 

A quantitative investigation of the time delay function 
is shown in Fig.[71 In the left panel, we plot At as function 
of the radial distance r for different values of the model 
parameters. Note that, as expected, the function At{r) has 
a maximum and two minimum thus meaning that three 
images are formed. In the right panels, r is held fixed and 
we plot At as function of /3 for different values of a. An 
image forming at r will be more and more delayed as a and 
(3 get higher. A detailed analysis of the dependence of At 
on the model parameters is quite complicate because of the 
many parameters involved (including the source position) 
so that will be not performed here. 

5.2. Lens ellipticity and contribution of the shear 

Real galaxies are hardly spherically symmetric. An easy 
way to take into account the ellipticity of the lens should 
be to insert it by hand in the lensing potential replacing 
the cylindrical dimcnsionless radius f with its elliptical 
counterpart rh — \/S^+QJcp with q the axial ratio. A 
still simpler approach is to add an external shear to the 
lensing potential since it has been shown that an elliptical 
potential ip(irii?) with an on axis shear 7s and axial ratio q 
produce the same image configuration as a pure elliptical 
potential with an axial ratio q' — g[(l— 7s)/(l-l-7s)]^^^ and 
no shear (|Witt 1996p . Moreover, an external shear could 
also mimic the effect of tidal perturbations from nearby 
galaxies or departures from a smooth distribution because 
of substructures in the primary lens. Motivated by these 
considerations, we add the shear term : 



ip^ 



1 



7sa^r cos2(6'- e'^J 



(53) 



to the lensing potential ip being (7s , 9^^ ) the shear strength 
and direction. 

The main effect of the shear is to alter the shape of 
the critical curves, defined as the loci in the lens plane 
where the magnification formally diverges. In particular, 
not considering the radial critical curve, the higher is 7s, 
the larger is the deformation of the tangential critical 
curve. In the case of no shear, this curve coincides with the 
Einstein ring, while, for 7s ^ 0, we have to distinguish two 
closed curve: the first one elongated along the direction of 
the shear itself (the critical curve) and the other one in 
perpendicular direction (the Einstein ring). Moreover, the 
larger is 7s, the lower will be the axial ratio h/a between 
the shortest and the longest axis of these ellipses. 
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Fig. 7. Left. Time delay function At (measured in days) as a function of the radial distance from the centre (measured 
in arcsec). Note that negative values of r corresponds to radial distances after a rotation of tt radians. Model parameters 
are set as (a,To,a;o) = (—0.3,5.0, 1.0) and three values of /3 are considered, namely (3 — 0.8 (solid), f3 = 1.0 (long 
dashed), f3 — 1.2 (short dashed). For the luminosity profile, we set 7 = 1, while the source position is rg — 0.3 arcsec. 
The extremals of the function corresponds to the images. Right. At as function of /3 for three values of a, namely 
a ~ —0.3 (solid), a = —0.5 (long dashed) and a = —0.7 (short dashed). The remaining model parameters are set as 
before, while we move the source at =0.1 arcsec and evaluate At at r = 1 arcsec. 
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Fig. 8. Image configurations for different values of the parameters : the dimension of the spot changes varying the 
parameters. Left. We set the source position = 0.1" and let /3 increasing (the higher is (3, the greater is the spot). 
Right. Here (3 is fixed at 1 and changes. The bold continue line is the critical curve and the corresponding caustic 
is the astroid (dashed bold curve at center); the internal lighter continue ellipse is the radial critical curve and the 
corresponding caustic is the dashed fighter ellipse. In both the panel we set a = —0.3, To — 5, 6s = 20°, 7^ = 0.15 
and 6*^. = 180°. 



5.3. Image formation 

As quoted above, the axisymmetric lens can produce only 
2 images with in addition a central (observationally not 
visible) one. When we distort the lensing isopotential con- 
tours by adding the shear term, the shape of tangential 
critical curves changes and the point - like caustic unfolds 
generating four cusps ( [Schneider, Ehlers fc Falco 1992^ . It 
is interesting to discuss the formation of the images as a 
function of some of the model parameters for the case with 
the shear term included. 

In Fig.m we plot the image configuration for some val- 
ues of the parameters and changing /3 or r^. We can see 
how the model can generate 3 or 5 images (2 or 4 plus 



a central image). In the left panel, we fix the source po- 
sition and change the values of f3 finding that the radial 
coordinate of the image increases with p. The right panel 
shows, for a fixed /?, how the number and the position of 
the images change with the source position r^. Increasing 
two of the images approach each other in direction of 
the tangential critical curve, until the source exits the re- 
gion delimited by the tangential caustic curve (the central 
astroid) when the two images merge and the system passes 
from a 5 - images to a 3 - images configuration. The other 
two curves in figure represents the radial critical curve and 
the corresponding caustic. 
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The number of images needs a brief remark. The odd- 
number theorem predicts that lenses with a smooth and 
not singular surface mass density k{R) which decreases 
faster than as |i?| — > necessarily forms an odd num- 
ber of images. For models with k diverging for i? — s- 0, 
we can use a generalized version of the previous theorem 
given by Evans & Wilkinson (1998). These authors con- 
sider models with an inner cusp and k cx (e.g., NFW 
models above described) and distinguish three different 
behaviours (see Ref. above for major details) 

1. if < n < 1 (weak density cusp), then we observe an 
odd number of images 

2. if 1 < n < 2 (strong density cusp), then we observe an 
even number of images 

3. if n = 1 (isothermal density cusp), then we can observe 
both an odd or even number of images 

Our model can predict an asymptotically decreasing (flat) 
K for 7 — 2 < a < 7— 1 (7— 1 < a < 7). When 
7 — 2 < a <7— lisn = — (a — 7-I-I) and it is sim- 
ple to verify that these models belong to the weak density 
cusp class, with < n < 1, thus an odd number of images 
always exist as observed in Fig. [5] for all the values of 7. In 
addition, when k is cored (for instance, if 7 — 1 < a < 7), 
then the original odd-number theorem holds and predicts 
also in this case an odd number of images. Thus, our model 
cannot reproduce strong cusp, predicting for all the pos- 
sible combinations of the lens parameters an odd number 
of images and the formation of both radial and tangential 
critical curves (the first one not observed in the case of 
strong cusp models) . Nearly all the observed lenses present 
two or four images, thus a discordance appears. However, 
observations and theoretical predictions can be reconciled 
if the central image is very close to the lens and highly de- 
magnified by the high central surface density, making this 
image difficult to detect. Actually, central images were ob- 
served in the 2-image lens PMN J 1632-0033 (Winn et al. 
2003, 2004), therefore the existence of these events can al- 
low to constrain the lens model and particularly the inner 
slope of the mass density. 

5.4. Comparison with commonly used lens models 

Having investigated in detail the lensing properties of our 
phenomenological model, it is worth comparing it with 
some of the most used models in literature. 

As a quite general and simple class of models, we con- 
sider the power law lensing potential : 

i;{R,0)^R''F{e,q,eg) (54) 

where < < 2 is the slope of the radial profile 
and F(9,q,9q) the angular shape function with {q,Oq) 
the axial ratio and the orientation of the isopoten- 
tial contours. Eq. ([54|) represents a generalization of the 
the pseudo - isothermal elliptical potentials (jKovner 1987| 
IKassiola fc Kovner 1995P that have been designed to 
be as similar as possible to the pseudo - isothermal 



mass distribution models (|Blandford fc Kochanek 1987|, 
IKassiola fc Kovner 1995^ . Given their easy manageability, 
they have been widely used in literature. For semiana- 
lytical and/or numerical treatments in modeling lensed 
quasars see, e.g., Cardone et al. (2001, 2002), Tortora et 
al. (2004). In addition, considering that k cx and 
following the previous Sec, power law with < v < 1 
{1 < < 2) have strong (weak) cusp. 

It is worth stressing that the models defined by Eq. ([5^ 
are characterized by a constant radial slope of the lens- 
ing potential. Obviously, this is not the case for our phe- 
nomenological model. Nevertheless, we can match the two 
models in the inner regions by setting a — 7 -I- 3 = i'. In 
particular, for 7 = 1 and C Oi our model reproduces 
the pseudo - isothermal trend for a G [—1,0] thus exclud- 
ing the region (0, 1], since for these values of a is > 2 
and the pseudo - isothermal models generate a non physi- 
cal surface mass density. It is possible to match the outer 
slope imposing a + P = v: also in this case the matching 
is not complete since when a + (3 > 1 our models are not 
defined, on the contrary the models in Eq. ((54|) holds also 
for V e (l,2j!l. 

6. Conclusions 

The rotation curves of spiral galaxies have been the first 
observational cornerstone of the dark matter theory. The 
nature and the distribution of this unseen component of 
galaxies are still obscure, but more and more evidences 
have been accumulated in favour of its ubiquitous pres- 
ence. Nonetheless, this picture is far to be free of its own 
problems since recent observations of line of sight veloc- 
ity profile of elliptical galaxies have furnished contrast- 
ing evidences favoring or disfavoring the presence of sig- 
nificant amounts of dark matter in the outer regions of 
these systems. In an attempt to face this problem we have 
employed here a phenomenological approach assuming a 
double power law expression for the global M/L ratio as 
T(7') = To(r/ro)"(l -I- r/ro)^ . This simple parametriza- 
tion reduces to constant M/L models for (a,/3) = (0,0), 
while gives an asymptotically increasing M/L ratio for 
a -|- /3 > as in models with dominating dark matter 
component. Adopting spherical symmetry and a versatile 
expression for the luminosity density j(r), it is possible 
to build a wide class of (effective) galaxy models start- 
ing from the mass profile M{r) = T(r)L(r). Varying the 
model parameters (a, [3) in the range determined by some 
physical considerations makes it possible to mimic a wide 
range of models having finite (a + /3 = 0) or formally 
infinite {a + (3 > 0) total mass, cuspy (a — 7 < 0) or 
cored (a — 7 = 0) density profiles, asymptotically flat 
(a + /3 — 1 = 0) or decreasing {a + f3 — 1 < 0) rotation 
curve, cuspy (7 — 2 < a < 7 — 1) or cored (7 — 1 < a < 7) 
surface mass density. 

* This circumstance is verified since we required that our 
model never has an increasing Vc (giving the constraint a + f3 < 
1), on the contrary models (|54p predict both decreasing and 
increasing Vc. 
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The main advantage of our parametrization of the 
M/L model is therefore the possibihty to reproduce a wide 
phenomenology resorting to a simple analytical expres- 
sion and relying on phenomenological considerations only 
thus avoiding any theoretical prejudice or systematic er- 
rors related to convergence problems in numerical N - body 
simulations of structure formation. As a preliminary and 
mandatory step, to allow a future comparison between 
theoretically predicted and observed properties of early 
type galaxies, we have here derived the main dynamical 
and lensing quantities of the model giving whenever pos- 
sible analytical expressions or illustrative plots. We have 
payed a particular attention to the asymptotic trends in 
the inner and outer regions since both these regimes are 
mostly probed by different kind of observations. On the 
one hand, measurement of the velocity dispersion profiles 
hardly extends to radii larger than the effective radius i?e 
and, indeed, only the central velocity dispersion is often 
at disposal. Such a mass tracer therefore mainly probes 
the inner region of the system so that knowing how our 
model fares for x « \ helps in comparing with this kind 
of data. In order to probe larger radii, one may resort to 
multiply imaged lensed quasars. First, the images typi- 
cally form at r ~ (1 — 2)i?e so that a larger region is 
probed. Moreover, the lensing potential and the fiux ra- 
tios also depend on the derivatives of the mass profile thus 
allowing to probe indirectly also outermost regions. 

We can individuate two main implications. Firstly, we 
have seen how the the full mass density profile and de- 
rived DM profile depend on 7 and therefore on the fitted 
luminous profile. Therefore, the presence of the baryons 
influences dark matter profile, since the parameter 7 that 
specifies the inner slope of the luminous profile enters in 
the expression of halo profile derived from our model. 
Secondly, we have compared our model with the more 
common generalized NFW profile: the asymptotic trends 
has been compared and we are able with suitable choices 
of a and /3 to recover the behaviour of gNFW model. The 
value of the inner trend of mass profile has generated a 
long debate in literature. The simulations seems to indi- 
cate an inner trend that scales like r^^ with 6 >1 (the so 
called cusped models) . On the contrary, many results from 
observational analysis, by rotation curves of low surface 
brightness, dwarf galaxies and spiral galaxies, dynamical 
studies of elliptical galaxies and gravitational lensing 
dMcGaugh fc de Blok 19981 [Binney fc Evans 2001 

IBorriello fc Salucci 20011 Ide Blok et al. 2001 



Treu fc Koopmans 2002b 



IBorriello et al. 2003P sug- 
gested the presence of internal core with 5 < 1 (cored 
models). This contradiction generates the so called 
cusp/core problem. Our model can reproduce both cuspy 
models with 5 >1 and cored ones with (5 < 1. In details, 
a < recovers the first models, i.e. the results from 
simulations and a > the seconde one, those from 
results of observations. Being able to reproduce all these 
different behaviours, our phenomenological model offers 
the perspective to homogenously analyze the observations 
and give an answer to the cusp/core problem. 



On a theoretical side, our approach could be further 
ameliorated. A key ingredient in our modeling is the lu- 
minosity density j{r) which we have assumed to be de- 
scribed by the class of Dehnen models. Although this 
gives a quite general expression able to reproduce well 
the observed surface brightness profile, the ideal route to 
follow should be to deproject an empirical function ob- 
tained by a direct fit to the photometric data. However, 
this approach should be repeated for each single galaxy 
and thus does not allow to investigate general properties. 
An intermediate possibility could be to use the depro- 
jection of the Sersic profile as input for the luminosity 
density and investigate the dynamical and lensing proper- 
ties of the corresponding model. Moreover, it has been 
shown that the dynamics of inner regions of elliptical 
galaxies is better described adding the contribution of the 
central supermassive black hole ( Baes fc Dcjong he 2004[ 
IMamon fc Lokas 2005al IMamon fc Lokas 2005b|l so that 
it is interesting to update our parametrization includ- 
ing this term into the mass profile. Finally, departures 
from spherical symmetry could also be investigated since 
real galaxies are likely to be moderately triaxial systems 
(|Franx et al. 1994p . 

Contrasting the model with the observations is, of 
course, the best strategy to both test the viability of our 
parametrization and constrain its parameters. Given the 
difficulties in obtaining radial velocity dispersion profile 
up to large radii, an ideal tool could be represented by 
lens galaxies. With more than 90 observed multiply im- 
aged quasars, gravitational lensing offers a unique dataset 
to put severe constraints on the slope parameters (a, (3) 
and hence suggest useful hints to solve the problem of dark 
matter in elliptical galaxies. Moreover, since lens galaxies 
span a range in redshift from low to intermediate, such an 
analysis also offers the possibility to investigate whether 
and how the characteristics of the dark matter content of 
early type galaxies evolve. 

A first qualitative result can be obtained considering 
the index of the surface mass density defined as indK(r) = 
dhiK,/ d\nr. In order to shape lensing events with their 
model independent analysis, Williams fc Saha (2000) have 
used as a constraint on the chosen lens model the condition 
indK(r) < —0.5 for values of r close to the images positions 
in strongly lensed quasars (see also ITortora et al. 2004p . 
Such a condition ensures that the image magnification is 
less than ~ 2 which is a quite plausible constraint well 
satisfied in some real lenses. Considering that typically 
images forms at (1 — 2)i?e, we have checked that the con- 
straints on indK(r) is satisfied by our model for all physi- 
cally reasonable choices of the slope parameters. Although 
this does not allow to further constrain (a,/?), this result 
suggests that our parametrization of the M/L ratio leads 
to a model that correctly reproduces this aspect of the 
observed lensing phenomenology. 

A detailed investigation of some interesting cases is 
outside our aim here and will be presented in a forth- 
coming paper. Here, we only note that, although the lu- 
minous component parameters (7, a) may be set from 
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the analysis of the surface brightness profile, severe de- 
generacies still exist among the slope parameters (a, (3) 
and the scaling quantities (Tq^xq)- As a consequence, 
it is likely that reliable constraints could be obtained 
only considering bright quadruply imaged systems so that 
the errors on the image positions is low and the num- 
ber of constraints higher. Moreover, it is preferable to 
limit the attention to isolated lens galaxies in order to 
avoid to introduce further unknown parameters to de- 
scribe the lens environment. Finally, it is highly desir- 
able to have a measurement of the central velocity dis- 
persion in order to break some of the degeneracies in- 
herent in lens modeling with the help of dynamical data 
dTreu fc Koopmans 2002al|Treu fc Koopmans 2"004[ ). 

As a final comment, we would like to stress again the 
need for a phenomenological approach to the problem of 
dark matter in elliptical galaxies. In order to overcome 
any theoretical prejudice, one should adopt a procedure 
that matches a large set of observational constraints with 
the minimal number of hypotheses. Adding gravitational 
lensing data to the dynamical one makes it possible to 
advance our knowledge on the observational side of the 
problem. Using phenomenological models like the one we 
have presented here could be the theoretical ingredient 
best suited to be coupled to the above data to analyze the 
question of dark matter in early-type galaxies. 

Acknowledgements. We warmly thank the anonymous referee 
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paper. 
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